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Abstract 


A single reversible solid oxide fuel cell (RSOFC) can perform dual functions: (1) as a solid oxide steam electrolyzer (SOSE) for hydrogen 
production and (2) as a solid oxide fuel cell (SOFC) for power generation. Thus, RSOFC can potentially offer a low-cost approach to support 
hydrogen economy. A modeling study has been conducted to analyze the important concentration overpotentials in both SOSE and SOFC modes 
of operation. The quantitative analyses show that in the SOSE mode, the hydrogen electrode is vulnerable to high concentration overpotential 
and limiting current density. Oppositely, in the SOFC mode, the oxygen electrode is vulnerable to above problems. If the SOSE and SOFC 
modes are considered separately, a RSOFC should be oxygen-electrode-supported and hydrogen-electrode-supported, respectively. For this reason, 
comprehensive analysis is very important to optimize the structure of the electrode-support to maximize the overall efficiency of a RSOFC 
performing dual functions. The modeling study signifies the difference between the SOSE and SOFC modes and provides insights in the operating 


mechanisms of RSOFC. The present model can be further extended to conduct more simulations for design optimization. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Hydrogen has been identified as an ideal energy carrier to 
support sustainable energy development. Hydrogen can be effec- 
tively produced by reforming fossil fuels, e.g. natural gas steam 
reforming [1—4]. However, the greenhouse gases produced will 
harm the environment and fossil fuels will be depleted beyond 
their economical usage in the foreseeable future. Development 
of renewable, clean, and economical hydrogen production tech- 
nologies to replace fossil fuel-based hydrogen production meth- 
ods is a key step towards a sustainable hydrogen economy [2]. 

Water electrolyzers integrated with photovoltaics or wind 
turbines will play an important role in renewable hydrogen pro- 
duction [5]. Among different types of electrolyzers, solid oxide 
steam electrolyzers (SOSE) are more advantageous due to their 
high efficiency [6-10]. Recently, the U.S. Department of Energy 
(DOE) initiated several projects on hydrogen production using 
SOSE [11-13]. 
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Since, a SOSE has the same electrodes and electrolyte as 
a solid oxide fuel cell (SOFC) but works in a reverse SOFC 
process, the concept of reversible solid oxide fuel cell (RSOFC) 
is feasible. That means a RSOFC can be operated in the SOSE 
mode to produce hydrogen fuel and the same RSOFC can later 
consume the hydrogen fuel to generate electricity in the SOFC 
operating mode [14—19]. The capability of dual functions makes 
RSOFC economically sound. 

In the present literature, some experimental studies on 
RSOFC are available but there is a lack of detailed modeling 
studies, which are indeed important for a better understanding 
of the operating mechanisms as well as for design optimization. 
Under the same operating condition (i.e., same temperature and 
same current density), a RSOFC working in either the SOSE or 
SOFC mode has identical ohmic and activation overpotentials. 
The overpotentials can be determined by the models reported 
previously [20-36]. On the other hand, the concentration over- 
potentials are different between the SOSE and SOFC modes 
because of the different gas transport mechanisms in the porous 
electrodes. Therefore, the concentration overpotentials are the 
sole factor responsible for the different RSOFC current-voltage 
(J-V) characteristics between the SOSE and SOFC modes. For 
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Fig. 1. Operating mechanisms of RSOFC in (a) SOSE mode; and (b) SOFC mode. 


this reason, the emphasis of this study is on modeling of the 
concentration overpotentials of a RSOFC. 


2. Model development 
2.1. Operating mechanisms 


The fundamental mechanisms of a RSOFC working in the 
SOSE mode are illustrated in Fig. 1a. Steam is fed to the porous 
hydrogen electrode (cathode). When an electrical potential, suf- 
ficient for water-splitting, is applied to the two electrodes, the 
water molecules diffused to the cathode-electrolyte interface are 
dissociated into hydrogen gas and oxygen ions. The hydrogen 
molecules produced are diffused out of the hydrogen electrode. 
The oxygen ions are transported through the electrolyte to the 
porous oxygen electrode (anode) and then oxidized to form oxy- 
gen molecules. Subsequently, oxygen, being the only substance 
in the mass transport in the porous oxygen electrode, is driven 
out by means of permeation. The net SOSE reaction can be 


expressed as 
H20 —> H; + 50>. (1) 


The working mechanisms of RSOFC in the SOFC mode with 
air as the source of oxidant are illustrated in Fig. 1b. Unlike the 
SOSE mode, oxygen and nitrogen, instead of pure oxygen, are 
present in the oxygen electrode. Thus, oxygen is transported by 
self-diffusion in the SOFC mode [23,24]. 


2.2. RSOFC cell potentials 


The J-V characteristics of RSOFC can be evaluated by the 
cell potentials related to the Nernst potential (E), activation over- 
potentials (jact,H, and Nact,o, ), ohmic overpotential (Hohmic), and 
concentration overpotentials (con, H) and Necon,o,) in the SOSE 
mode, 


SOSE SOSE 
Vsose = E+ Nact, H + Nact, O2 + Nohmic + Neon „H2 F Ncon,O2> 


(2) 
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and in the SOFC mode, 


SOFC _ „SOFC 
Ncon, H2 Ncon, O2 : 


(3) 


The Nernst potential can be determined by the Nernst equa- 
tion. The activation overpotentials can be determined by 
the Butler-Volmer equation [17,18]. The ohmic overpotential 
dependent of the current density, electrolyte thickness, and oper- 
ating temperature can be obtained by [20] 


Vsorc = E — Nact, Ha — Nact,0. — Nohmic 


10300 
Nohmic = 2.99 x 1075JL exp (F) (4) 


where T (K) is absolute temperature and L (pm) is the thickness 
of the electrolyte. 


2.3. Concentration overpotentials in the SOSE mode 


The concentration overpotential is caused by the resistance 
to the transport of reactant species approaching the reaction 
sites and the transport of product species leaving the reaction 
sites. In the SOSE mode, the concentration overpotentials can be 
expressed in terms of the gas concentration difference between 
the electrode surface and the electrode-electrolyte interface as 
[23,24,26] 


I 0 

„SOSE = RA i Ph, Pho (5) 

con, H2 2F PO, Pl o 
and 

1/2 
RT Pl 

SOSE O 

Ncon,O. = OF In (3 ) (6) 
2 


where P° and P! are the partial pressures at the electrode surface 
and electrode—electrolyte interface, respectively, and subscripts 
H20, H2, and O2 represent the steam reactant, hydrogen and 
oxygen products, respectively. 

It is assumed that the electrochemical reactions take place at 
the electrode—electrolyte interfaces [23,24]. Each participating 
component undergoes a steady-state mass transport in the porous 
electrode [26,40,41], 


VN;=0 (7) 


where N; is the rate of mass transfer of species i. In the SOSE 
porous hydrogen electrode (cathode), the mass transport mech- 
anism is based on diffusion as the gas molar generation rate 
is equal to gas molar consumption rate. Both Fick’s model 
and Dusty Gas Model (DGM) are frequently used to describe 
gas transport within porous media. DGM is based on Stefan 
Maxwell equations and includes Knudsen diffusion. However, 
no analytical expression can be obtained and numerical method 
must be employed. In comparison, Fick’s model is the sim- 
plest model that can describe the gas transport characteristics 
effectively. Furthermore, analytical solutions can be obtained by 
using Fick’s model. In the present study, the diffusion process, 


which is driven by concentration gradient or partial pressure 
gradient, is determined by Fick’s model [40], 

N, = DI BOP) DE ACP) 

RT ox RT ox 

where P is the operating pressure, x is the depth measured from 
the electrode surface, Pj, yj, and pe are the partial pressure, 
molar fraction, and effective diffusion coefficient of species i, 
respectively. Substituting Eq. (8) into Eq. (7), the equation gov- 
erning the mass transport of species i can be obtained, 


a 
Ox 


ER oa (8) 


Ds 3(P;) 
RT ox 


=0 (9) 


where i represent H20 and H3 at the cathode. It is noted that for 
advanced porous electrodes of which reactions can take place 
in the porous media, the micromodel developed by Chan and 
co-workers can be employed [37-39]. 

At the cathode-electrolyte interface (x = dp,, the thickness 
of cathode), the diffusion rate of H20 towards the interface is 
equal to the H20 consumption rate under steady state condi- 
tion. As H20 consumption is governed by the electrical current 
density J, the mass transfer rate can be expressed as 


J 
N = —. 10 
mo = 55 (10) 


Substituting Eq. (10) into Eq. (8), the boundary condition for 
the governing Eq. (9) becomes 


PHO eS (11) 
Ox x=duyy 2F Dio 


The Dirichlet boundary condition can be applied for the cathode 
surface, 


Piol x= = Pho: (12) 
Solving Eq. (9) with the above two boundary conditions, the 
partial pressure of H20 at the cathode-electrolyte interface can 
be obtained, 


Pho = Pho — AE L ii (13) 
Defi, UF 


As the steam molar consumption rate is equal to the hydrogen 
molar generation rate, the cathode total pressure (P = Py, + 
Pp,o) is constant. Therefore, 


I I 0 RT J 
Py, = P— Pho =P- Pho = Det zp 
2 
RT J 


= Ph, + sd. (14) 
7 Dibo 2F 


Substituting Eqs. (13) and (14) into Eq. (5), the concentration 
overpotential of the hydrogen electrode can be obtained, 


sage. RT es Fee + (JRTdy, /2 FE) 


n H 7 
con, H2 2F Pi, (Pho =- (JRTdp, /2FD}f 0)) 
_ RT P 1+ (JRTdy, /2F Dis PG) (15) 
2F \ 1- (JRTdn,/2FDf oP) 
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The diffusion in the porous SOSE electrodes is mainly based 
on two mechanisms, namely, molecular diffusion and Knudsen 
diffusion. Molecular diffusion is the dominant mechanism if the 
pore size is much larger than the mean free path of the molecular 
species. In this case, the molecule—molecule interaction governs 
the diffusion process. On the other hand, if the pore size is much 
smaller than the mean free path of the species, the molecule—pore 
wall interaction dominates over the molecule—molecule inter- 
action. Thus, Knudsen diffusion becomes an important mech- 
anism. In most porous structures, both mechanisms are sig- 
nificant. The effective diffusion coefficient of steam can be 
expressed by combining these two diffusion mechanisms using 
the Bosanquet formula [26,40,41], 


1 == ( 1 4 1 ) (16) 
Dio € \Dyo-Hw Dmo,k 


where &/e is the ratio of cathode tortuosity to porosity, 
&/ (e Dp,0-np,) is the reciprocal of the effective molecular diffu- 
sion coefficient for a HyO-Hp binary system, and &/(¢Dy,0,x) 
is the reciprocal of the effective Knudsen diffusion coefficient 
for steam. 

For the Knudsen diffusion, as gas molecules frequently col- 
lide with the walls of the pores, the transport of molecules can 
be modeled using the kinetic theory [26,40,41], 


4 SRT 
Dp,o,k = 37 Ta (17) 
2 


where r is the mean pore radius and My, is the molar weight 
of H20 (18 g mol~!). The molecular binary diffusion coefficient 
Dy,0-H, can be obtained from the Chapman-Enskog theory of 
ideal gas [41], 


1 1 
Dy,0-H, = 0.00135 ( 7 + Mao 
2 2 


1/2 T3/2 
) (18) 


2 
POOH 2p 


where My, is the molar weight of H2 (2 g mol7!); OH0,H> is the 
mean characteristic length of species H20 and Ho, and 2p is a 
dimensionless diffusion collision integral. The analytical values 
of 2p and oy,0,H, can be obtained from references [26,41]. 

At the oxygen electrode (anode), O2 is the only gas 
present in the porous electrode layer. The O2 generated at the 
electrolyte—anode interface causes a pressure gradient in the 
porous electrode. Therefore, O2 is transported in the electrode 
by permeation, instead of diffusion. Darcy’s law is the most 
widely used model to characterize gas viscous flow in porous 
media, depending on the gas dynamic viscosity and the porous 
microstructure. Employing Darcy’s law, the oxygen flux can be 
determined [22,42], 


(Po/80,) — (Po/80,) — P8,) exp |(RT/4F)(J80,do,/ DË po)| 


where jz is the dynamic viscosity of O2, which can be deter- 
mined by sixth-order polynomial functions developed by Todd 
and Young [43]; and Bg is the flow permeability, which can be 
determined by the Kozeny—Carman relationship as [44,45] 


3 
E 
B; = —— —30ry. 20 
g A ) (20) 
Substituting Eq. (19) into Eq. (7), the governing equation for 
oxygen transport in the porous anode can be obtained as 


0 Po, Bg dPo, -o0 
ox RT ox | ` 


(21) 


At the anode-electrolyte interface (x = do,, the thickness 
of anode), the rate of O2 transporting away from the interface, 
which equals to the O2 generation rate under steady state con- 
dition, is also governed by the electrical current density, 


No. = -75 (22) 


The boundary condition of Eq. (21) for the anode can be obtained 
as 


Po, - = . (23) 


x=do, 


The Dirichlet boundary condition can be applied for the anode 
surface, 


Pos |0 = Po» (24) 


Solving Eq. (21) with boundary conditions expressed in Eqs. 
(23) and (24), the partial pressure of O2 at the anode—electrolyte 
interface can be obtained, 


2 JRTudo, 
PL, = y (PL + =. (25) 


2FBy 


The anode concentration overpotential can thus, be expressed as 


0 \2 
sose _ RT In Ve) + (JRTudo, /2FBg) 
"con,O2 = 4F Pe 
2 


(26) 


2.4. Concentration overpotential in the SOFC mode 


The following modeling equations developed by Chan and 
co-workers [23,24] can be employed to determine the concen- 
tration overpotentials of a SOFC, 


(1 — (RT/2F))(Jdy, / Dé P8, ) 
(1 + (RT/2F))( Jd, / Dif, P§,0) 


SOFC RT 


Ncon, H2 ia oF n 


| » 7) 


RT 
SOFC 
"con,O2 = AF In po > (28) 
02 
and 
p£ 
Po, Bg 60, = —— k _ 29 
No, =——"_2V Po, 19 © = (29) 
02 RT  ” oe DÈ x + D-n 
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where Po, DĚ, DÆ p and DË n are the pressure of 
the oxygen electrode, effective oxygen diffusion coefficient, 
effective Knudsen diffusion coefficient of oxygen, and effec- 
tive oxygen-nitrogen binary diffusion coefficient, respectively. 
These diffusion coefficients can be determined by the same 
method described in previous section. 


2.5. Computation procedures 


The J-V relationship of RSOFC operating in the SOSE mode 
is implicitly formulated by Eqs. (2), (4), (15), and (26). The 
analytical Butler-Volmer equations from references [17,18] are 
directly employed to model activation overpotentials. Using 
J as the independent variable, the numerical values of corre- 
sponding ohmic overpotential, activation overpotentials, con- 
centration overpotentials, as well as the cell potential can be 
computed by direct substitution of values into the right-hand 
side of the equations. For the SOFC mode, the J-V relationship 
can be obtained by the concentration overpotentials expressed 
by Eqs. (27) and (28) and the cell potential expressed by Eq. 
(3). 


3. Results and discussion 
3.1. Model validation 


The above model has been derived to determine the J—V char- 
acteristics of RSOFC performing dual functions. The model for 
the SOFC mode developed by Chan and co-workers [23,24] 
is well accepted (Eqs. (27) through (29)); therefore, no vali- 
dation is needed. On the contrary, the newly developed model 
for the SOSE mode should be validated by comparison with 
experimental data. Momma et al.’s experimental work in SOSE 
hydrogen production has been selected because the laboratory 
setup and test procedures are clearly reported in the literature 
and the details can facilitate the theoretical simulation [46]. In 
Momma et al.’s experiments, the J-V characteristics of hydro- 
gen production by electrolyte-supported planar SOSE discs were 
measured. The electrolyte, cathode, and anode are made of 
yttria stabilized zirconia (YSZ), nickel- YSZ (Ni- YSZ), and lan- 
thanum strontium manganite (LSM), respectively. The thickness 
of the electrolyte, cathode, and anode are 1000 wm, 100 um, 
and 100 um, respectively. The tests were conducted under a 
constant pressure of | atm, steam molar fractions between 20 
and 60%, and various operating temperatures from 1173 to 
1273 K. In the theoretical simulation, the values of the input 
parameters are summarized in Table 1 [23,24,26]. The effects 
of operating temperature and current density on cell poten- 
tial are shown in Fig. 2. Fig. 3 shows the dependence of 
the cell potential on the steam molar fraction for J equal to 
2000 Am~? and T equal to 1273 K. The cell potential decreases 
as the steam molar fraction increases. From Figs. 2 and 3, 
the good agreement between the theoretical simulation results 
and the experimental data by Momma et al. [46] serves as 
a thorough validation of the theoretical model for the SOSE 
mode. 


Table 1 

Values of input parameters in the present study 

Parameter Value 
Operating temperature, T (K) 1073 
Operating pressure, P (atm) 1.0 
Exchange current density at the hydrogen electrode (A/m?) 5300 
Exchange current density at the oxygen electrode (A/m?) 2000 
Electrode porosity, € 0.3 
Electrode tortuosity, Ẹ 6.0 
Average electrode pore radius, r (um) 0.5 
Electrolyte thickness, L (um) 1000 
Electrode thickness, d (um) 

For model validation 100 
For parametric analyses 50, 500 


3.2. Concentration overpotentials of hydrogen electrode 


Based on a parametric modeling analysis, the concentration 
overpotentials of the RSOFC hydrogen electrode operated in 
both SOSE and SOFC modes are compared in Fig. 4. The results 
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Fig. 2. Comparison of simulation results with experimental data [46] for model 


validation—effect of operating temperature and current density on RSOFC poten- 
tial in SOSE mode. 
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Fig. 3. Comparison of simulation results with experimental data [46] for model 
validation-effect of steam molar fraction on RSOFC potential in SOSE mode. 
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Fig. 4. Comparison of SOSE mode with SOFC mode: concentration overpoten- 
tials of hydrogen electrodes. 


generally show that the concentration overpotential decreases 
with increasing reactant gas molar fraction, i.e. H20 for SOSE 
and H2 for SOFC. This finding is consistent with the previous 
studies that higher cell potential can be achieved by increasing 
hydrogen molar ratio at the inlet of a SOFC [26]. It is also found 
that the RSOFC hydrogen electrode has higher concentration 
overpotential in the SOSE mode than in the SOFC mode. The 
limiting current density in the SOSE mode is generally lower 
than that in the SOFC mode. These phenomena can be explained 
by the diffusion of reactant gases in the porous hydrogen elec- 
trode. As H20 has much higher molecular weight than H2, the 
Knudsen diffusion coefficient of H2O is smaller than that of 
H2. Consequently, the RSOFC has a smaller effective diffusion 
coefficient in the SOSE mode than in the SOFC mode (see Eqs. 
(16) and (17)). The SOSE concentration overpotential of the 
hydrogen electrode becomes a more crucial limiting factor. 


3.3. Concentration overpotentials of oxygen electrode 
Fig. 5 compares the oxygen-electrode concentration overpo- 
tentials between the SOSE and SOFC modes. It is found that the 
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Fig. 5. Concentration overpotentials of oxygen electrodes of RSOFC: SOSE 
mode vs. SOFC mode. 


tentials is insignificant at low current density below 1400 A/m?. 
As the current density increases, the SOSE concentration over- 
potential increases moderately while the SOFC concentration 
overpotential increases considerably resulting in a limiting cur- 
rent density. 

The SOFC limiting current density occurs because the dif- 
fusion rate of reactant O2 in the porous oxygen electrode 
is reaching the limit as the O2 concentration at the oxygen- 
electrode—electrolyte interface is approaching zero at high cur- 
rent density [23,24]. The low diffusion coefficient due to high 
molecular weights of O2 and N2 yields low effective diffusion 
coefficient and thus, a low limiting current density. On the other 
hand, there is no limiting density current for the oxygen elec- 
trode in the SOSE mode of operation because the transport of 
oxygen ions (reactant) through the ion-conducting electrolyte is 
not limited by the porous electrode. 


3.4. Electrode-supported RSOFC 


The concentration overpotentials of SOFC have been exten- 
sively investigated. It is found that hydrogen-electrode supports, 
generally more advantageous than oxygen-electrode supports, 
yield low overpotential loss and high power output [23,24]. This 
finding is certainly applicable to RSOFC operating in the SOFC 
mode. 

The concentration overpotentials of RSOFC operating in the 
SOSE mode are studied here in detail with the present theo- 
retical model. In practice, a cell consists of a thick layer of 
electrode (typically 500 um) that can provide structural support, 
while the other electrode is a thin layer (50 um) to minimize 
overpotential loss. The characteristics of different electrodes 
(thick or thin, hydrogen- or oxygen-) are shown in Fig. 6. A 
thick hydrogen electrode suffers seriously from high concen- 
tration overpotential and low limiting current density. For the 
oxygen electrode, a thick layer only causes a slight increase in 
the concentration overpotential. More importantly, the oxygen 
electrode in the SOSE mode does not suffer from limiting cur- 
rent density. The results imply that an oxygen electrode-support 
performs better than a hydrogen electrode-support in the SOSE 
mode. 
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Fig. 6. Concentration overpotentials of RSOFC electrodes in SOSE mode. 
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Based on the above findings, the selection of an electrode 
support can greatly affect the overall performance of a RSOFC. 
A hydrogen-electrode support is favorable to the SOFC mode 
while an oxygen-electrode support is favorable to the SOSE 
mode. Therefore, the details of both SOSE and SOFC operating 
conditions should be carefully considered in design of RSOFC 
on a case-by-case basis. For example, if a RSOFC is specified 
mainly to generate high electricity, ahydrogen-electrode support 
should be chosen in the preliminary design. On the contrary, if 
a RSOFC is mainly used for hydrogen generation, an oxygen- 
electrode support is recommended. 


4. Conclusion 


A theoretical model has been successfully developed to inves- 
tigate the concentration overpotentials of a RSOFC, which can 
perform dual functions as a SOSE for hydrogen production and 
as a SOFC for power generation using hydrogen fuel. The analyt- 
ical results signify the difference between the SOSE and SOFC 
operating modes due to different gas transport mechanisms. The 
theoretical model provides insights of the RSOFC operations. 
The model can also be effectively applied to subsequent opti- 
mization of the RSOFC design for maximum overall energy 
efficiency. 
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